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ABSTRACT 

Cosmological applications of the “redshift - angular size” test require knowledge of the linear size of 
the ’’standard rod” used. In this paper, we study the properties of a large sample of 140 milliarcsecond 
compact radio sources with flux densities measured at 6 cm and 20 cm, compiled by Gurvits et 
al.(1999). Using the best-fitted cosmological parameters given by Planck/WMAP9 observations, we 
investigate the characteristic length l m as well as its dependence on the source luminosity L and 
redshift l m = IL^( 1 + z) n . For the full sample, measurements of the angular size 9 provide a tight 
constraint on the linear size parameters. We find that cosmological evolution of the linear size is small 
(|n| ~ 10 -2 ) and consistent with previous analysis. However, a substantial evolution of linear sizes 
with luminosity is still required (/? ~ 0.17). Furthermore, similar analysis done on sub-samples defined 
by different source optical counterparts and different redshift ranges, seems to support the scheme of 
treating radio galaxies and quasars with distinct strategies. Finally, a cosmological-model-independent 
method is discussed to probe the properties of angular size of milliarcsecond radio quasars. Using the 
corrected redshift - angular size relation for quasar sample, we obtained a value of the matter density 
parameter, Q rn = 0.292to;ggp, in the spatially flat ACDM cosmology. 

Subject headings: quasars: general galaxies: active - radio continuum: galaxies - cosmology 


1. INTRODUCTION 

The redshift - angular size data have provided 
a us eful method to pr o be c o smological parame¬ 
ters (IcTu erra. Dalv fc Wa r] 120001 |Vishwa Wmal 120011 : 


iLima fc Alcanizl 120021 IChen fc Ratra 12003^ 7 since this 


relation is directly related to the angular diameter dis¬ 
tance. Powerful radio sources constitute a population 
which can be observed up to very high redshifts, reach¬ 
ing beyond feasible limits of supernova studies. Over a 
few past decades considerable advances have been made 
to investigate the redshift - angular size relation in ra¬ 
dio sources for the pu rpose of cosmological studies, in- 


1995; 

Buchalter et al. 1998c Gurvits, Kellerman & Frey 

1999: 

Guerra. Dalv & Wan 200(1 Zhu & Fuiimoto 

2002: 

Podariu et al. 2003; Jackson 2004j; Barai & Wiita 

2006, 


been measured for different t ypes of radio sources, suc h 
as extended FRIIb g alaxies dDalv fe Diorgovskil [200.'} ). 
radio lo ud quasars dBuchalter et all 1199811 and radio 
galaxies dCuerra. Dalv fc Waiill200(il h In a similar spirit, 
by using radio observations of the Sunyaev-Zeldovich ef¬ 
fect (SZE) together wi t h X-r a y emission of galaxy clus ¬ 
ters. iDe Filippis et al.l (|2005l) : iBonamente et al.1 (|2006d 
extensively explored the angular diameter distances at 
different redshifts. 

In fact, in order to break inherent degeneracies be¬ 
tween cosmological parameters every alternative method 
of restricting these parameters is desired in mod¬ 
ern cosmology. Consequently, there have been nu¬ 
merous atte mpts to use compact radio sources for 
this purpose fVishwakarma 200 l| ILima fe; Alcanizl 120021 : 
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IZhu &; Fuiimotol l2002t IChen fe Ratral [2003h . In these 
studies the analysis was carried out on 12 binned data- 
points with 12-13 compact sources per bin. One of 
the major uncertainties was the typical value of the 
linear size l. In order to obtain cosmological con¬ 
straints, so me authors chose to fix l at certain spe¬ 
cific values (iVishwakarnial |2001tJLima ii &Alcanid 12002 : 
IZhu fc Fuiimotoll2002t l . while lChen fc Ratral (|2003l i chose 
to include some range of values for l and then marginal¬ 
ized over it. 

It is obvious that cosmological application of the “red¬ 
shift - angular size” data requires good knowledge of 
the linear size of the “standard rod” used. The pos¬ 
sibility that source’s linear size depends on the source 
luminosity and redshift should be kept in mind. In 
particular it still remains controversial whether com¬ 
pact radio sources are indeed “ true” standard rod s 
(jCurvits. Kellerman fc Frevl 119991 : iVishwakarmal I2001H . 
Applying th e po pul ar parametr ization l m = IL^{\ + 
z ) n , I Gurvits. Kellerman fe Frevl (j!999T) claimed that 
by excluding sources with extreme spectral indices 
and low luminosities their compiled data has been 
minimized for the confounding by these two effect s 
(jCurvits. Kellerman fc Frevl 119991 : IVishwakarmal l20C)Tlf . 
However, their results were obtained under assump¬ 
tion of a homogeneous, isotropic universe without cos¬ 
mological constant (Ha = 0). On the other hand, 
from a large body of recent as tronomical observations , 
such as Union2 SNe la dataset (jAmanullah et a,lJl20 1 01 1 . 
the CMB observation from the W ilkinson Micro wave 
Anisotropy Probe (WMAP9) (jHinshaw et al.ll20T3i ). and 
the BAO distance ratios from the spectroscopic Sloan 
Digital Sky Survey (S P SS) d ata r elease7 (DR7) galaxy 
sample dPadmanabhan et al.l I2012T) , no convincing evi¬ 
dence for deviations from the concordance ACDM model 
has been established. More recently, Planck, the third- 
generation space mission following COBE and WMAP, 
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has recently released its first cosmological results based 
on measurements of the CMB temp eratu re and lensing- 
potential power spectra (lAde et al.l 120141) . All of them 
strongly indicate the existence of an exotic component 
called dark energy, which represents more than 70% of 
the total energy of the universe and serving as a driv¬ 
ing force of the cosmic acceleration. Latest investiga¬ 
tions specifically to study the prope r ties of dark energy 
were c a rried out bvlYu fe Zhul (1 20 1 11) : I C a.o. Liang fe Zhul 
Eouh: IC^ fe Zhul (120121): ICao. Covone fe Zhul (120121): 


iCao et al.l (l2012l)~ Cao. Zhu fc Zha ol (12012 1): iPan et al.1 


(1201211: ICa o i^Liangl ()20131) : il.iao. Pan feZhul "<|20T3^ : 

ICao fc Zhul ( 2014D . Having this in mind, properties 
of compact radio sources should be readdressed with 
present angular size data and taking into account a reli¬ 
able cosmology based on current precise observations. 

In this paper, we will reconsider issues associated with 
angular sizes of radio sources under the assumption of 
ACDM cosmological model. Specifically, we will study 
the characteristic length, the “angular size - luminos¬ 
ity” and “angular size - redshift” relations for the com¬ 
pact structure in quasars and radio galaxies assuming 
ACDM cosmological model as the background, which 
is well supported by observations. In Section 2 and 
3 we briefly describe our sample, its construction and 
methodology of subsequent analysis. Results with the 
full sample and several sub-samples are presented in Sec¬ 
tion 4. Cosmological-model-independent constraints on 
the compact source parameters and their cosmological 
application are discussed in Section 5. Finally, we sum¬ 
marize the conclusions in Section 6. 


2. OBSERVATIONAL DATA 

Our goal is to better constrain the parameters model¬ 
ing compact radio sources, i.e. their linear size scale l , 
along with luminosity and redshift dependence of their 
metric linear length l m . By “better” we mean ob¬ 
tained using the best currently available cosmological 
model. To this end, we have considered the angular 
size d ata fo r milli arcsecon d radio sources compiled by 
iGurvits. Kellerman fe Frevl (119991) . This data set was a 
larger sample of source s than used bv lKellermannl (I1993D 
or by Wilkin son et al.1 (119981) and w i th more com plete 
structural data than used bv IGurvitsI (| 19931 119941) . 

All 330 sources included in this comprehensive com¬ 
pilation were imaged with very-long-baseline interfer¬ 
ometry (VLBI) at 5 GHz with resolution of ca. 1.5 
mas. They included: 1) 79 compact radio sources as¬ 
sociated with acti ve galaxies and quasars considered in 
iKellerman n (1993); 2) s ources described by the Caltech 
- Jodrell Bank group dXu et, al.l 119951: iHenstock et, al.1 
Il995t iTavlor et al.l 11994 I1990D : 3) sources discussed 
in other works as well as observat i ons of high-redshif t 
{z > 3) quasars dFrev et al.1119971 iParaei et al.l I1998D . 
IGurvits. Kellerman fe Frevl (I1999D then reduced this 
original compilation of 330 sources down to 145 sources, 
with spectral index (—0.38 < a < 0.18) and to¬ 
tal luminosity ( Lh 2 > 10 26 WHz -1 ). They claimed 
that the former criterion helps excluding compact 
sources with inverted spectrum and relatively large 
steep spectrum, while the latter tends to minimize 
the possible dependence of linear size on luminosity 
(IGurvits. Kellerman fe Frevlll999l) . 

Full information about all the 145 sources that remain 


afte r the aforementioned selecti on, can be found in Table 

1 of IGurvits. Kellerman fc Frevl (|1999D . including source 
coordinates, redshifts, optical counterpart, angular size, 
spectral inde x, and tot al flux densities at 6 cm or 20 cm. 
Following iKellermannl ( 1993 ). the characteristic angular 
size of each source was defined as the distance between 
the strongest component (referred to as the core) and 
the most distant component which had the peak bright¬ 
ness greater than or equal to 2% of the peak brightness 
of the core. We emphasize here, that in order to imple¬ 
ment multi-frequency “0-z” tests in our analysis (see the 
next section for details), we further restricted the final 
sample to 140 data points with measured flux densities 
at both 6 cm and 20 cm. The final sample, which covers 
the redshift range 0.031 < 2 < 3.89 and does not show 
jet-like structure for any system, contains a wide class 
of extragalactic objects including 112 sources identified 
as quasars, 18 radio galaxies, and 10 BL Lac objects 
(blazars). 

Radio galaxies we use in this work are located within 
the redshift range of 0 < z < 0.9. This means that, 
using the “redshift - angular size” test, we are able to 
constrain properties of active galaxies and their evo¬ 
lution up to 2 ~ 0.9. Main motivation of study¬ 
ing the milliarcsecond radio structures in q uasars stems 
from their potential usefulness in cosmo logy (|Kellermannl 
119931 : IGurvits. Kellerman fe Frevl 1 1 999D . Moreover, the 
large number of quasars is also beneficial for studying 
structural properties of milliarcsecond radio structures 
at high redshifts. Finally, inclusion of 10 sources with BL 
Lac objects as counterparts allows us to show that their 
structural properties are similar to the known quasars at 

2 > 3. 

However, the “angular size - redshift” test requires sta¬ 
tistically complete and well-characterized (homogeneous) 
sample. Because our list includes sources corresponding 
to different optical counterparts at different redshift as 
described above, so besides the full combined sample we 
will also consider separately seven sub-samples. Three of 
them are defined by optical counterparts: radio galaxies, 
quasars and blazars. Next subsamples are defined by re¬ 
striction to four redshift ranges: 2 < 0.5, 0.5 < 2 < 1.0, 
1.0 < 2 < 2.0 and 2 > 2.0. 


3. METHOD 

According to lSandagel ( 1988 ). the angular size-redshift 
relation for a rod of intrinsic length l m can be written as 


e{z) 


D a (z ) 


( 1 ) 


where l m is the metric linear size, Da is the angular di¬ 
ameter distance at redshift 2 . Following the_ phenomeno¬ 
logical model first proposed in IGu rvitsI (Il994 ) and later 
discussed in IGurvits. Kellerman fe Frevl (|1999D , the pro¬ 
jected linear size of a source is related to its luminosity 
L and redshift 2 as 

lm=K^) P {l+z) n ( 2 ) 

where l is the linear size scaling factor representing the 
apparent distribution of radio brightness from the peak 
down to the level of its 2% in the available sample of 
VLBI images. It is a parameter defined by the practical 
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Table 1 

Summary of constraints on the metric linear size parameters obtained with the full sample and six sub-samples (see text for definitions). 
Source of cosmological priors and the wavelength of measured flux are given in brackets. 


Sample (Cosmology+Flux density) 

l (pc) 

P 

n 

Full sample (Planck+Se) 

Full sample (WMAP9+Se) 

Full sample (Planck+S^o) 

l = 25.42 ± 3.62 
l = 24.91 ±3.59 
l = 27.07 ±3.95 

p = 0.169 ±0.025 

P = 0.170 ±0.026 
p = 0.183 ±0.027 

n = -0.021 ±0.139 
n = -0.009 ±0.141 
n= -0.090 ±0.145 

Radio galaxy (Planck+S'e) 

Quasar (Planck+Se) 

BL Lac (Planck+Se) 

L — ^y-OO_20.05 
l = 25.96 + 4.14 
l = 20.73jlii;7g 

P = 0.242 ± 0.064 
p = 0.203 ± 0.034 

P = 0.139 ±0.107 

n = 0.142 ±0.670 

n = -0.023 ±0.153 
n = -0.315 ±0.895 

Sub-sample (z < 0.5) (Planck+Se) 
Sub-sample (0.5 < z < 1.0) (Planck+Se) 
Sub-sample (1.0 < z < 2.0) (Planck+Se) 
Sub-sample (z > 2.0) (Planck+Se) 

l = 44.60 + 21.80 

/ — 21 fi e i+ 21 - 20 

l — zi.uo_ 13 20 
l = 63.85tSi;J! 

1 = 10.49±! 4 3 f 

p = 0.284 ± 0.076 

P = 0.009 ± 0.091 
p = 0.215 ±0.083 
p = 0.128 ±0.187 

n = 1.228 ± 1.072 
n = -0.335 ± 1.115 

n = -0.880 ±0.600 

n = 0.518 ±0.638 


limitation of dynamic range of VLBI data, i.e., a higher 
sensitivity of VLBI observations would enable estimates 
of angular sizes down to lower values of brightness, re¬ 
sulting in turn in a differ ent value of the characteristic 
linear size (jGurvi tHHH). l 0 is the normalizing lumi¬ 
nosity taken to be equal to I0 28 WHz -1 in our analysis. 
Parameters /? and n represent the dependence of the lin¬ 
ear size on source luminosity and redshift, respectively. 
The first parameter - /3 is related to physics of a com¬ 
pact radio emitting regions, while the parameter n mim¬ 
ics three physical effects: (1) cosmological evolution of 
the linear size with redshift; (2) dependence of the lin¬ 
ear size on the emitted frequency; and (3) an impact of 
sources bro adening d u e to scatter ing in_the propagation 
medium (jGurvits. Kellerman fc Frevi 11999 !. The latter 
effect is not important for our sample with the lowest 
emitted frequency of 5 GHz (corresponding to z = 0). 
The distinction between the former two effects require 
multi-frequency 9 — z tests. 

The luminosity L of radio sources is estimated from 
their measured flux density S 0 b s • So the radio luminosity, 
assuming isotropic emission, reads: 


L = S ° hA \^f (3) 

(l±z) 1+ “ v ' 

where Sobs is the observed flux density, Dl is the lu¬ 
minosity distance, and a is the spectral index (Sobs ^ 
v a ). In general, for sources at cosmological distances, 
k-correction must be applied to the spectral index a of 
the source. The angular diameter distance D A and the 
luminosity distance Dl at redshift z are related to each 
other through the so-called distance duality relation 

D l /D a (1 + z)~ 2 = 1 (4) 

which is a fundamental relation in observational cosm ol¬ 
ogy and initiated a lo t of studies, e.g. dCao fe Liangl 
l201lt IGao fcZh^l20nl) . 

The above equations imply that, if we could have a reli¬ 
able knowledge of cosmological model parameters which 
therefore allow to calculate D A or Dl at different red- 
shifts, then we would get stringent constraints on the 
range of parameters, l, /?, and n describing compact ra¬ 
dio sources. Theoretical expression for the angular di¬ 
ameter distance D A {z) (expressed in Mpc and assuming 


flat FRW metric) reads 


D a {z', Urn) 


3000 / 1- 1 r z dz' 

(1 + 2 ) Jo E(z'' : Q m ) 


( 5 ) 


where h is the dimensionless Hubble constant, E(z; f2 m ) 
is the dimensionless expansion rate, which - in the case 
of flat ACDM model - depends on redshift z and matter 
density parameter fl m in the following way: 

E 2 (z\ Cl m ) = Q m ( 1 + z) 3 + (1 — fi m ). (6) 

For the purpose of our analysis, theoretical D A (z) has 
been calculated by using the best-fit matter density pa¬ 
rameter given by Planck Colla borati on: fi m _ = 0.315 ± 
0.017 and h = 0.673 ± 0.012 (lAde et al.l[20Tl . Even 
though Planck results are the latest ones, we also in¬ 
clude the data from the Wilkinson Microwave Anisotropy 
Probe 9 year data release (W MAP9), i.e. = 0 .279 ± 
0.025 and h = 0.700 ±0.022 (jHinshaw et all 20131. The 
h = 0.700 value is also used the cosmological application 
of the cosmological-model-independent method is dis¬ 
cussed in Section 5. In order to determine the parameters 
of compact radio structures, we preformed Monte Carlo 
simulations of the posterior likelihood C ~ exp(—y 2 /2) 
using routines available in CosmoMC package. As a 
prior, we assumed a conservative 20% Gaussian uncer¬ 
tainty of the observed angular size. 


4. RESULTS AND DISCUSSIONS 

In this paper, we focused our attention on the con¬ 
straints on the parameters (l, /3, and n) characterizing 
compact radio sources obtained from different samples, 
i.e. the full N = 140 sample as well as several sub¬ 
samples determined from different selection criteria. The 
results are summarized in Table 1. 


4.1. Estimates on the full sample 

As we already remarked, measured flux density at dif¬ 
ferent bands could bear the information about physical 
conditions in active galactic nuclei - a feature common 
to all types of sources we used. Performing fits on the 
data comprising flux at 6 cm, we obtained the following 
best-fit values and corresponding ler (more precisely 68% 
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confidence level) uncertainties 

l = 25.42 ±3.62 pc, 
/3 = 0.169 ±0.025, 
n = -0.021 ±0.139. 


Then, using the flux densities at 20 cm we obtain the 
following best fit 


l = 27.07 ±3.95 pc, 
/3 = 0.183 ±0.027, 
n = -0.090 ±0.145. 


Marginalized lcr and 2<r contours of each parameter ob¬ 
tained at different bands are shown in Fig. [T] It is clear 
that the parameter degeneracies are consistent with each 
other, as can be seen in the contours obtained from flux 
densities measured at 6 cm and 20 cm, respectively. Best 
fitted values obtained for different wavelength are differ¬ 
ent, but they agree within lcr. 

It is obvious that, for well resolved compact sources, 
measurements of 9 provide tighter estimates of the lin¬ 
ear size parameters (l, /3, n). More importantly, our full 
sample analysis has also yielded improved constraints on 
the meaningful physical parameters: /? and n. We found 
that best-fitted value of the parameter n is a small num¬ 
ber: \n\ ~ 10~ 2 , slightly negative, but 68% Cl contains 
zero in any case i.e. our results are consistent with no 
evolution of l m with z. This suggests that, contrary to 
the case of extended radio sources, central engine pow¬ 
ering compact radio sources is likely to be controlled 
by a limited number of physical parameters (mass of 
central black hole, accretion rate) and may therefore 
be less subject to evolutionary effects. On the other 
hand, we found that, for 140 sources satisfying lumi¬ 
nosity selection criterion Lh? > 10 2e WHz~ x , substan¬ 
tial evolution of linear size with luminosity is still re¬ 
quired. Comp ared w ith previous resu lts obtained on the 
same data (iGurvits. Kellerman fc Frevlll999l) , our results 
show that, improved, more rigorous quantitative analy¬ 
sis supports the existence of “linear size - luminosity” 
relati on. T he conc lusion tha t 0 + n > —0.15 given by 
IGurvits. Kellerman fc Frevl (1 9991 ) does not contradict to 
the findings of the present work. The best-fit parame¬ 
ters of the phenomenological dependence Eq.© under 
the modern cosmological model are different from those 
obtained with the_ c lassi cal Einste in - de Sitter model 
used bv IGurvits. Kellerman fc Frevl (1 999 ). The values 
of the two best-fit parameters of the phenomenological 
formula obtained here, namely, /? and n , if confirmed by 
future ’’angular size-redshift” studies, would offer addi¬ 
tional constraints for cosmological tests based on angular 
sizes of extragalactic sources. 

As we stressed above, the assumption of currently best 
available cosmological model ACDM was the source of 
improvement concerning estimates of l, /3, and n. Hence 
their values depend on the cosmological parameters used. 
Therefore, besides assuming flat ACDM model with pa¬ 
rameters coming from Planck observations, we also con¬ 
sidered WMAP9 results for comparison. In this case, the 



Figure 1. Constraints on compact source parameters obtained 
from the full sample, based on different flux density measurements 
at 6 cm and 20 cm, respectively. 

best fit is 


1 = 24.91 ±3.59 pc, 

/3 = 0.170 ±0.026 
n = -0.009 ±0.141. 

Marginalized probability distributions for each param¬ 
eter and marginalized 2D 68% confidence contours are 
presented in Fig. [2] Comparing constraints based on 
Planck and WMAP9 observations, we see that confidence 
regions of l, j3 , and n are almost the same, hence our re¬ 
sults and discussions presented above are robust. We 
remark here that, considering the WMAP9 and Planck 
data are consistent with the accuracy sufficient to the 
comparison with the “0-z” test, it is not surprising the 
regression results of the “0-z” test in combination with 
WMAP and Planck are compatible in the framework of 
ACDM cosmology. 

4.2. Estimates on sub-samples 

In Table 1 and Fig. [3J131 we show the results of fit¬ 
ting three parameters, l, f3 and n on seven sub-samples 
described in Section [2] 

We note that the ranges of l and f} parameters for 
quasars (l = 25.96 ± 4.14 pc, /3 = 0.203 ± 0.034) 
are marginally close to estimates obtained for compact 
structures in BL Lac objects (l = 20.73l 2 'j ! '7g PC P = 
0.139 ± 0.107). Rather weak dependence of the char¬ 
acteristic size on redshift, i.e. the range of the param¬ 
eter n for quasars ( n = —0.023 ± 0.153) is in agree¬ 
ment with the estimate obtained for BL Lac sources 
(n = —0.315 ± 0.895) within lcr. On the other hand, 
luminosity dependence (/3 = 0.242 ± 0.064) and weak 
redshift dependence (n = 0.142 ±0.670) are both present 
in radio galaxies. The best-fit values of /3 and n for 
this sub-population are significantly different from the 
corresponding quantities of quasars or BL Lac sources. 
Consequently, our results imply the existence of physical 
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Figure 2. Constraints on compact source parameters obtained 
from the full sample, assuming cosmological parameters inferred 
from Planck and WMAP9 data. 


Figure 4. Constraints on compact source parameters obtained 
from sub-samples defined by four different redshift bins. 



Figure 3. Constraints on compact source parameters obtained 
from three sub-samples with different optical counterparts. 

differences between galaxies and quasars at the milliarc- 
second scale. To some extent, this conclusion supports 
the scheme of treating radio galaxies and quasars with 
distinct strategies. We must keep in mind that similarity 
or difference in (/3, n) parameters for radio sources with 
different types of optical counterparts, might reveal sim¬ 
ilar or different physical processes governing the radio 
emission of compact structures. 

This tendency could also be found in fits performed 
on four sub-samples with different redshift bins. From 
Fig. 2] we find that: (1) Constraints on all the param¬ 
eters coming from the low-redshit sub-sample (z < 0.5) 
are substantially different from those obtained with other 


Figure 5. Constraints from the N = 26 quasar sub-sample, with 
observational luminosity distance derived from Union2 SN la data. 


sub-samples. This can be explained by the fact that that 
low redshift sub-sample is dominated by radio galaxies. 
(2) For the sub-sample with redshift range 0.5 < z < 1.0, 
the “no-evolution” model (0 = n = 0) is still included 
within lcr confidence regions in (/3, n) parameter plane, 
whereas a substantial evolution of linear sizes with lumi¬ 
nosity is still required for the other three sub-samples. 

As we remarked above, sub-samples defined by red¬ 
shift ranges are confounded with types of optical counter¬ 
parts. Therefore a stratified analysis taking into account 
redshift ranges and the source type would be desirable. 
However, our sample is too small to achieve this. 
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Figure 6. Cosmological fits on the ACDM model form the cor¬ 
rected angular size - redshift relation for quasars (Blue line) and 
joint data sets combined with BAO observations (Red line). 

5. COSMOLOGICAL-MODEL-INDEPENDENT 
CONSTRAINTS ON COMPACT SOURCE PARAMETERS 

In the previous section, we discussed the constraints 
on the model parameters of compact structure in ra¬ 
dio sources using theoretical expression for angular di¬ 
ameter distances and assuming best currently available 
ACDM cosmology. In this section, we propose cosmolog¬ 
ical model independent approach. Namely we derive the 
angular diameter distance for our radio sources from ob¬ 
served luminosity distance of S N la in the Union2 compi¬ 
lation (lAmanullah et al1l2010l) . This provides us a nat¬ 
ural way to calibrate the properties of angular size of 
milliarcsecond radio sources. 

In order to place cosmological-model-independent con¬ 
straints on /, (3, and n , one should first perform pair¬ 
wise matching of radio sources and SN la almost at the 
same redshift. Since the sample size of the SN la is much 
larger than that of the radio sources, we bin the observed 
Da (inferred from the Union2 data points) according to 
the following criterion: A z = \z ra dio ~ zsn < 0.005 
(|Cao fc Liangi [2011 ). As a result we obtain a sample 
of 42 observational angular diameter distances Da de¬ 
rived from the supernova data covering the redshift range 
0.031 < z < 1.34. However, not all of them could be 
used in the cosmological-model-independent method. As 
we already discussed in Section 2 and Section 4, full 
N = 140 sample is not statistically complete and ho¬ 
mogeneous. Moreover, sources corresponding to differ¬ 
ent optical counterparts may have distinct “angular size 
- redshift” relation. Therefore, in cosmology indepen¬ 
dent analysis we limited ourselves to radio sources with 
quasars as counterparts. There are two reasons sup¬ 
porting this choice. First, our sample is dominated by 
quasars. Second, as we have seen, cosmological evolution 
of the linear size for quasars is very small, so we can as¬ 
sume n = 0. Therefore, we finally used SN la to derive 
observational angular diameter distance for 26 quasars. 

Fitting results of the compact source linear size param¬ 


eters (l, fi) are shown in Fig. [5} with the best fit 

l = 26.70 ±2.91 pc, 

P = 0.158 ±0.040. 

One can see that the results derived from the 
cosmological-model-independent analysis agree very well 
with the best-fitted parameters determined from theo¬ 
retical cosmological distances for the TV = 112 quasar 
sample. 

Having performed cosmological-model-independent 
analysis, we can consider cosmological implications of 
the corrected redshift - angular size relation. Using the 
best-fitted l and /3 parameters (with their la uncertain¬ 
ties) obtained from the model-independent analysis to 
the full quasar data and performing “angular size - red¬ 
shift” test assuming flat ACDM model, we are able to 
get the observational constraint on the matter density 
parameter. The result is f \ m = 0.292+ O '.o90 and its pos¬ 
terior probability density function is shown in Fig. [6] We 
see that it is in agrement with the value obtained from 
Planck observations within la range around the central 
value. Moreover, our analysis result is fully compati¬ 
ble with that obtained from the previous study of pecu¬ 
liar velocities of galaxies, fl m = 0.30)to’o7> which is the 
only alternative metho d sens itive exclusively to matter 
density (Feldman et al. 2003). Based on the 12 binned 
data-points with 12-13 compact sources per bin, it was 
found that Friedmann mod el with a vanishing A is not 
the best-fit cosmology (IVishwakarmalf^OO l). This result 
supported the necessity to include dark energy i n the cos¬ 
mological model. Then [Lima fc Alcanizl ( 2002 ) used the 
same binned data to place constraints on the flat XCDM 
cosmology (including dark energy with constant equation 
of state w) with fixed physical length l m ~ 20 — 30 h~ l 
for the radio sources. They demonstrated that the flat 
ACDM model with = 0.20 and l m = 22.6Ti -1 is the 
best fit to these milliarcsecond radio source data. The 
potential of using the same sample to study other cos¬ 
mological models including dark energy with constant 
or t i me-variab le equation of state was also discussed in 
iChen fe Ratral (120031 ). More recently, by applying an as- 
tropliysical model to quantify the behavior of compact 
radio sources as standard rods and considering possible 
selection effects, I.Tacksoni (120041 ) gave the best-fit param¬ 
eter = 0.24)to p 7 for the flat ACDM model from the 
original data set (jGurvi ti l 1994 ). We find that the con¬ 
straints resulting from our analysis are consistent with 
the previous works. However, because we used the cur¬ 
rently favored cosmological model and performed a cos¬ 
mological model independent check, our results could be 
useful as hints for priors on l, /3 and n parameters in 
future cosmological studies using compact radio sources. 

It has been known for some time that cosmologi¬ 
cal parameters can be more stringently determi ned us¬ 
ing additional and complementary data (iCao fe Zhul 
I2014D . Therefore, we combine the compact radio sources 
“angular size - redshift” test with the latest data on 
Baryon Acoustic Oscillations (BAO). More specifically 
we add the BAO data from the Sloan Digital Sky 
Survey (SP SS) data r e lease 7 (DR7) corresponding to 
z = 0.35 dPadmanabhan et all 120121 ). SDSS-III Baryon 
Oscillation Spectrosco pic Survey (BOSS) at z = 0.57 
(lAnderson et all I2012D , the clustering of WiggleZ sur- 
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vey rtBlake et al.ll2012f ) at z = 0.44,0 60,0. 73 and 6dFGS 
survey at z = 0.10 (|Beutler et al.1 l201llb Likelihood 
distribution function for the fl m parameter in ACDM 
model constrained by BAO and the compact structure 
in quasars is also plotted in Fig. [6] Obviously, fl m is 
more tightly constrained with joint data set, with the 
best-fit parameter f 1 m = 0.308 ± 0.019. This also agrees 
with that obtained from Planck observations very well. 
Moreover, we find that cosmological constraint with the 
112 quasars is well consistent with the joint statistical 
analysis. 

6. CONCLUSION AND DISCUSSION 

In this paper, we explored the properties of a sample of 
140 milliarcsecond compact radio sources with measured 
angular sizes. Metric linear size of compact sources is 
usually parameterized as l m = IL@( 1 + z) n . Using the 
best available cosmological model parameters given by 
the Planck/WMAP9 observations, we investigated the 
elements of l m — its characteristic length l as well its 
dependence on the source luminosity L and redshift z. 
In the full sample, we found that measurements of 9 pro¬ 
vide tighter estimate of the source linear size parameters. 
Small cosmological evolution of the linear size |n| ~ 10~ 2 
is consistent with previous analyses, while a substantial 
evolution of linear sizes with luminosity is still required 
(/3 ~ 0.17). H owever, t h e conclusi on that 0 + n> —0.15 
given bv lGurvits. Kellerman fc Frevl (Il999f) does not con¬ 
tradict to the findings of our work. 

Furthermore, by dividing the full sample into seven dif¬ 
ferent sub-samples given the source redslrifts and their 
optical counterparts, we obtain the following results: (1) 
The range of the parameters /3 and n for quasars are 
close to the the estimates obtained for compact struc¬ 
tures in BL Lac sources within lcr. (2) Both luminosity 
dependence and weak redshift dependence are present in 
radio galaxies. The best-fit values of /3 and n for this 
sub-population are significantly different from the corre¬ 
sponding quantities of quasars or BL Lac sources. (3) 
Closeness or difference of parameter values for different 
types of counterparts, might reveal the similar or dif¬ 
ferent physical processes governing the radio emission of 
compact structures. This tendency could be also found 
from the constraints obtained with the sub-samples lo¬ 
cated at different redshift bins: (1) Constraints on the 
parameters with the low-redshit sub-sample (z < 0.5) 
are essentially different from those obtained with other 
sub-samples. (2) For the sub-sample with redshift range 
0.5 < 2 < 1.0, the “no-evolution” model (/3 = n = 0) is 
still included at lcr confidence region in the (3 — n plane, 
whereas a substantial evolution of linear sizes with lumi¬ 
nosity is still required for other three sub-samples. 

Finally, we studied the properties of angular size of 
milliarcsecond radio quasars with a cosmological-model- 
independent method, and then we derived the constraints 
(from the corrected quasar sample) on the spatially flat 
ACDM cosmology. The obtained value of matter density 
parameter, f 1 m = 0.292^°’^, agrees very well with the 
previous results obtained on the same “0 - z” sample 
and other recent astrophysical measurements including 
Planck observations. 

Therefore, our analysis indicates that, the radio source 
size seems to be dependent on the source luminosity, i.e. 
the sources are not “true” standard rod. This is incon¬ 


sistent with their model previously discussed in the lit¬ 
erature, in which this dependency has been minimized 
by discarding low values of luminosities and extreme val¬ 
ues of spectral indices. However, in order to differentiate 
observational selection effect from intrinsic luminosity- 
dependence, we still need multi-frequency VLBI obser¬ 
vations of more compact radio sources with higher sen¬ 
sitivity and angular resolution. 

As a final remark, we point out that the sample dis¬ 
cussed in this paper is based on VLBI images observed 
with various antennas configurations and techniques for 
image reconstruction. Our analysis potentially suffers 
from this systematic bias and taking it fully into account 
will be included in our future work. Moreover, the sta¬ 
tistical results are obtained with VLBI images observed 
at frequency of 5 GHz. Since the parameter n may re¬ 
flect possible dependence of the linear size on the emitted 
frequency, multi-frequency “9 - z” tests should also be in¬ 
cluded in the future work. 
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